First Measurement of and z/g Events in an Off-Axis Horn-Focused Neutrino Beam 
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We report the first observation of off-axis neutrino interactions in the MiniBooNE detector from 
the NuMI beamline at Fermilab. The MiniBooNE detector is located 745 m from the NuMI pro- 
duction target, at 110 mrad angle (6.3°) with respect to the NuMI beam axis. Samples of charged 
current quasi-elastic ffj, and Ve interactions are analyzed and found to be in agreement with expecta- 
tion. This provides a direct verification of the expected pion and kaon contributions to the neutrino 
flux and validates the modeling of the NuMI off-axis beam. 

PACS numbers; 25.30.Pt, 13.15. +g, 14.60.Lm, 14.60.Pq 



Conventional neutrino beams from high-energy proton 
accelerators serve as important tools for studying neu- 
trino characteristics and the fundamental properties of 
matter involving interactions of neutrinos. Such beams 
typically arise from the two-body decays of n and K 
mesons produced by a proton beam impinging upon a 
nuclear target. The mesons leave the target with a sig- 
nificant angular divergence. The flux of neutrinos in such 
a wide band beam at distance d from the meson decay 
point and at an angle 6 with respect to the parent meson 



direction is given by 

where 7 is the Lorentz boost factor of mesons To 
obtain a more intense neutrino flux, it is essential to fo- 
cus the mesons produced in the target. To accomplish 
this, neutrino experiments such as MiniBooNE 2] and 
MINOS 13] use focusing magnetic horns to direct the 
mesons toward downstream detectors. The energy of j/^s 
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FIG. 1: Plan and elevation views of the NuMI beamline with 
respect to the MiniBooNE detector. The MiniBooNE detec- 
tor is located 745 m from the NuMI production target, at 
110 mrad (6.3°) with respect to the NuMI beam axis. The 
length of the NuMI target hall is 45 m; the length of the 
evacuated decay volume is 675 m, and the distance from the 
target to the Minos Near Detector is 1040 m. 



from two-body decays is given by 



1 - 



Ett.K 



1 + 7^ tan^ 



(2) 



where mTr,K {Et^^k) is the mass (energy) of the tt, K 
parent, and is the muon mass. Brookhaven experi- 
ment E889 proposed Q an off-axis beam because, at a 
suitable off-axis angle 0, the neutrino flux is confined to 
a relatively narrow band of energies, which is useful in 
limiting backgrounds in searches for the oscillation tran- 
sition v^. Future neutrino oscillation searches by 



experiments plan to use off- 



the T2K and NOi/A 
axis horn-focused beams. 

The MiniBooNE detector, located at an angle of 
110 mrad (6.3°) with respect to the NuMI beam axis 
(see Fig. [T]), provides a unique opportunity to perform 
the first measurement of neutrino interactions from an 
off-axis horn-focused beam. In addition to demonstrating 
the off- axis beam concept, the measurement verifies the 
predicted fluxes from tt/K parents in the NuMI beam, 
and probes the off-axis intrinsic v^, contamination, re- 
quired for future — > appearance searches. 

The NuMI beam points toward the MINOS Far De- 
tector, located in the Soudan Laboratory in Minnesota. 
Neutrinos are produced by 120 GeV protons incident on 
a carbon target. In the period studied here, the beam 
intensity was up to 3 x 10^^ protons on target per spill 
at a typical repetition rate of 0.48 Hz. Positive tt and 
K mesons produced in the target are focused down the 
decay pipe using two magnetic horns. Neutrinos from 
two-body decays of pions are more forward directed than 
those from kaons due to the difference in rest mass of 
the decaying mesons. As a result, the off-axis component 
coming from pions is suppressed relative to the kaon com- 
ponent. Decay in flight of poorly focused pions can only 
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FIG. 2: Comparison of the predicted NuMI off-axis and on- 
axis fluxes including all neutrino species. The off-axis flux is 
separated into contributions from charged tt and K parents. 



occur close to the NuMI target, since they are stopped by 
shielding around the target area. The NuMI beam also 
provides a large sample of Vf, events in the MiniBooNE 
detector. The v^s result primarily from the three-body 
decay of kaons and thus have a wider range of energies. 
The stability of the neutrino beam is monitored using the 
muon monitors at the end of the decay pipe and the MI- 
NOS Near Detector. The direction of the neutrino beam, 
its intensity, and its energy spectrum were all found to 
be very stable over the data taking period [1] . 

Detailed GEANT3-based ^ Monte Carlo (MC) sim- 
ulations of the beam, including secondary particle pro- 
duction, particle focusing, and transport, are performed 
to calculate the flux as a function of neutrino flavor and 
energy. The yield of pions and kaons from the NuMI tar- 
get is calculated using the FLUKA cascade model 0. 
The beam modeling includes downstream interactions 
in material other than the target that produce hadrons 
decaying to neutrinos. These interactions are modeled 
using a GEANT3 simulation, configured to use either 
GFLUKA or GCALOR cascade models. The 
NuMI neutrino flux at the MiniBooNE detector is shown 
in Fig. O Pions (kaons) produce neutrinos with average 
energies of about 0.25 GeV (2 GeV). 

These neutrinos are detected in the MiniBooNE detec- 
tor [l3| which is a 12.2 m spherical tank filled with 800 
tons of pure mineral oil. The detector triggers on every 
NuMI beam spill and the detector activity is recorded in 
a 19.2 window beginning about 1 jjis before the start 
of a ~10 fjis wide spill. The stability of the NuMI beam 
was confirmed by studying the rate of neutrino events 
in the MiniBooNE detector as a function of time. The 
time and charge of photomultiplier-tubes (PMT) in the 
detector are used to reconstruct the interaction point, 
event time, energies, and particle tracks resulting from 
neutrino interactions. Neutrino interactions in the de- 
tector are simulated with the NUANCE event generator 
package with modifications to the quasi-elastic (QE) 
cross-section as described in [12[. Particles generated by 
NUANCE are propagated through the detector, using a 
GEANT3-based simulation which describes the emission 
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of optical and near-UV photons via Cherenkov radiation 
and scintillation. Neutrino induced events are identified 
by requiring the event to occur during the NuMI beam 
spill, after rejection of cosmic ray muons and muon de- 
cay electrons For the sample satisfying these selec- 
tion criteria, NuMI neutrinos are predicted to interact 
via charged-current (CC) QE scattering (39%), CC single 
pion production (31%), neutral current (NC) single pion 
production (14%), multi-pion production (9%), deep in- 
elastic scattering (4%), and other interactions (3%). The 
predicted event composition is v^^ : v^j^ : '■ Ve 0.81 : 
0.13 : 0.05 : 0.01. 

The data set analyzed here corresponds to 1.42 x 10^° 
protons delivered to the NuMI target from June 22, 2005 
to March 2, 2007. The MC in all cases has been normal- 
ized to this number of protons. There is a 2% uncertainty 
in the number of protons on target. Neutrino interac- 
tions are identified with the likelihood-based algorithm 
used in 

The high rate and simple topology of CCQE events 
provides a useful sample for understanding the j/^ spec- 
trum and verifying the MC prediction for Vf, production. 
The identification of CCQE events is based upon the 
detection of the primary stopping muon and the associ- 
ated decay electron as two distinct time-related clusters 
of PMT hits, called 'subevents': Vfi+n — > n~+p, n~ — > 
+ Ufj^ + De- We require the first subevent to have a 
reconstructed position within 5 m of the center of the 
detector. The decay electron requirement substantially 
reduces CC single 7r+ contamination because most CC 
TT"*" events contain a second decay electron from the 7r+ 
to /Lt decay chain. Additional rejection of non-v^ CCQE 
events in the sample is achieved by a requirement on the 
relative /i and e likelihoods, maximized under a given par- 
ticle hypothesis, log(Le/L^) < 0.02. This selection cri- 
terion is 24% efficient in selecting j/^ CCQE candidates, 
resulting in a 69% pure CCQE sample. The most sig- 
nificant background contribution to the i^^ CCQE sample 
results from CC single 7r+ production (78%) where the 
7r+ is undetected. A total of 17659 data events pass this 
ly^j. CCQE selection criteria, compared to the prediction 
of 18247±3189 in the 0.2 < < 3.0 GeV range; the un- 
certainty includes systematic errors associated with the 
neutrino flux, neutrino cross-sections, and detector mod- 
eling. 

The flux uncertainties include particle production in 
the NuMI target, modeling of the downstream inter- 
actions, and kaons stopped in the NuMI beam dump. 
The flux uncertainty also includes the uncertainty aris- 
ing from possible misalignment of the target, the focusing 
horns and the shielding blocks although this was found to 
have a small effect on the off-axis neutrino flux 1^, 14 1. 
The tt/K yields off the target were tuned to match the 
observed neutrino event rates in the MINOS Near De- 
tector where the same meson decays produce signif- 
icantly higher energy neutrinos. Such tuning has a neg- 



TABLE I: Observed and predicted //^-like events in two en- 
ergy bins. The predicted events are further separated into 
contributions from kaon and pion parents of neutrinos. 
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10734 ±104 


11169 ±1989 
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3534 


0.9-3.0 


6925 ± 83 


7078 ±1354 


1884 


5194 



ligible effect on the off-axis beam at MiniBooNE. How- 
ever, the difference between the tuned and the untuned 
tt/K yields was treated as an additional systematic ef- 
fect. Further details of systematic uncertainties consid- 
ered may be found in though it is important to note 
that the magnitudes of these systematic uncertainties due 
to the flux are smaller in the off-axis case. The cross- 
section uncertainties are quantified by varying the un- 
derlying model parameters constrained by either exter- 
nal or Booster neutrino beam (BNB) data. Uncertainties 
in the parameters describing the optical properties of the 
MiniBooNE detector are constrained by external mea- 
surements of the oil properties and by fits to calibration 
samples of events in the BNB data sample '5] . 

Reconstructed I'fj. CCQE event kinematics include the 
muon energy, i?^, and muon angle with respect to the 
neutrino beam direction, 0^. For both the data and MC, 
0^ is approximated assuming that all neutrinos arise from 
meson decays at the NuMI production target. In reality, 
mesons decay along the NuMI beamline so that the aver- 
age decay distance from the target is about 70 m. How- 
ever, given the geometry of the beamline with respect to 
the detector (see Fig. [T]), such an off-axis angle change is 
well within the angular resolution of the detector 2°). 
Based on these reconstructed quantities, the neutrino en- 
ergy. El,, is calculated assuming two-body kinematics 



2MpEf, 



Mr, 



{El 



mV) cos 6^ 



(3) 



where Mp is the proton mass. The Ei, resolution of NuMI 
neutrinos in MiniBooNE is ~ 12% at 1 GeV. The E„ dis- 
tribution of selected i^^ CCQE events is shown in Fig. [3l 
along with the MC prediction, separated into contribu- 
tions from target pions, target kaons, and non-target 
sources. About 50% of the events in the CCQE sample 
originate from parents produced in non-target materials. 
Predicted pion and kaon contributions, in two energy re- 
gions, are given in Table HI Systematic uncertainties of 
the predicted event rates are given in TablelHl The agree- 
ment between data and the prediction of the neutrino flux 
from tt/K parents indicates that the NuMI beam mod- 
eling provides a good description of the observed off-axis 



v^, flux in MiniBooNE. 



A measurement of the target kaon flux is performed 
using the CCQE sample. For this purpose the MC 
prediction is divided into four templates: neutrinos from 
kaons produced in the NuMI target (shown in Fig. [3]), 
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FIG. 3: Reconstructed E^, distribution of CCQE events. 
The band indicates the total systematic uncertainty associ- 
ated with the MC prediction. The prediction is separated 
into contributions from kaon parents produced in the NuMI 
target, pion parents produced in the NuMI target, and kaon 
and pion parents produced in non-target materials. 



FIG. 4: Mass distribution of NC n candidates for data 
(points) and MC (solid histogram). The dashed histogram is 
the subset of predicted events with at least one true vr". Pre- 
dicted non-vr'' backgrounds are either from and (dash- 
dotted line) or and i>e (dotted line) interactions. Kaon 
parents contribute 84% of the events in this sample. 



TABLE II: Systematic uncertainties of the predicted event 
rate in the full energy range and in two bins, for CCQE 
i/fi and i/e samples. 



Event 


Systematic Error 




E^[GcV] 




Sample 


Component 


0.2-3.0 


0.2-0.9 


0.9-3.0 




Flux [%] 


6.9 


7.2 


9.0 




Cross-sec. [%] 


15.7 


15.9 


16.2 




Detector [%] 


3.2 


3.6 


4.6 




Total [%] 


17.5 


17.8 


19.1 




Flux [%] 


9.8 


8.5 


11.9 




Cross-sec. [%] 


14.6 


14.2 


15.6 




Detector [%] 


8.5 


10.0 


8.9 




Total [%] 


19.5 


19.5 


21.8 



neutrinos from pions produced in the NuMI target (also 
shown in Fig. [3]) , neutrinos from kaons produced in non- 
target materials (i.e. downstream), and neutrinos from 
pions produced in non-target materials. Neutrinos from 
target kaons dominate the event rate in the Ei, > 1.2 
GeV region. The predicted target kaon fraction in this 
energy range is 83% of the selected sample. In order to 
get a clean measurement of the target kaon contribution, 
events with neutrino energies above 1.2 GeV are selected 
and a MC template fit to the data is performed. In the 
fit, the predicted target kaon component includes uncer- 
tainties associated with neutrino cross-section model and 
detector modeling; other MC components include these 
sources of error as well as flux uncertainties. After fit- 
ting the target kaon component of the MC prediction to 
the data, with the other MC components fixed to their 
initial values, the fit yields a kaon flux that is 1.14 +/- 
0.22 compared to the initial MC prediction. The fit has 
a x^/NDF = 9.8/10. Therefore, the measured flux of 
kaons from the NuMI target is consistent with the predic- 
tion derived from the FLUKA cascade model [1]. This 
measurement is a valuable as a check of kaon production 
which can provide a background in searches for j/^ — > 



transition in the NuMI beam. 

The CCQE events consist of a single subevent of 
PMT hits (ve+n -> e^+p). In 8% of j/^ CCQE events the 
/i~ is captured on carbon, resulting in a single subevent. 
These events are removed with an energy dependent re- 
quirement on the likelihood ratio, log{Le/ Lfj_). The ma- 
jority of the remaining background is NC tt*^ events with 
only a single reconstructed electromagnetic track that 
mimics a i^e CCQE event. To test our NC tt" predic- 
tion, a clean sample of NC tt*^ events is reconstructed, as 
shown in Fig. |4l This sample demonstrates good agree- 
ment between data and MC. About 30% of the events 
in the NC tt*^ sample originate from parents produced in 
non-target materials. The majority of tt*' events in the 
Ve CCQE sample are rejected by requirements on the re- 
constructed tt" mass and the electron to pion likelihood 
ratio, applied as a function of visible energy. 

A source of low energy events arises from the decay 
of stopped kaons in the beam stop at the end of the 
NuMI beamline, which is under the MiniBooNE detector 
(see Fig. [T]), 83 m from its center. Given the kinematics 
of stopped kaon decay, all VeS from this source will have 
visible energies {Ee) below 200 MeV. A requirement Ee > 
200 MeV effectively removes this source. A total of 780 
data events pass all of the CCQE selection criteria. 
The MC prediction is 660±129 with a i^e CCQE efficiency 
of 32% and purity of 70%. About 38% of the events in the 

CCQE sample originate from parents produced in non- 
target materials. The corresponding energy distribution 
is shown in Fig.[5l and the uncertainties on the predicted 
event rate are given in Table |TT1 To facilitate further 
comparison, the low and high energy regions are divided 
at 0.9 GeV, and the numbers of data and MC events in 
these two regions provided in Table IIIII The data with 
El, < 0.9 GeV are systematically above the prediction at 
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FIG. 5: Reconstructed E^, distribution of Ve CCQE can- 
didates. The prediction is separated into contributions from 
neutrino parents. The band indicates the total systematic un- 
certainty associated with the MC prediction. Kaon parents 
contribute 93% of the events in this sample. 



TABLE III: Observed and predicted i^e-like events in two 
energy bins. The predicted events are further separated into 
intrinsic i^e (and i/e) and (and i?^) components. 
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0.2-0.9 


496±22 


401±78 


311 


90 


0.9-3.0 


284±17 


259±56 


231 


28 



the 1.2 (J level. 

The results of the measurements described here show 
that reliable predictions for an off-axis beam can be 
made. However, it should be noted that the MiniBooNE 
experiment was not constructed as an off-axis detector 
for NuMI, but rather an on-axis detector for the BNB. 
The uncertainties in the neutrino flux presented here are 
substantially higher than might be expected in a long- 
baseline off-axis neutrino experiment, due to the close 
proximity of the MiniBooNE detector to the decay pipe 
and beam dump and the hadronic interactions therein 
which produce neutrinos. Consideration of such items 
will be important for future experiments being proposed 
or executed [3] • 



In summary, we have presented the first observation 
and analysis of neutrino interactions with an off-axis 
horn-focused neutrino beam. The agreement between 
data and prediction in the and Ve CCQE samples 
demonstrates good understanding of both pion and kaon 
contributions to the beam. This represents a successful 
demonstration of an off-axis neutrino beam at 110 mrad 
and provides a clear proof-of-principle of the off-axis 
beam concept planned for use in future neutrino experi- 
ments. 
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